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Governing Equations and Solution Methodology 

In the present work, a generalized form of the thin-layer Navier-Stokes equations is used to model the 
turbulent, viscous flow encountered in a synthetic jet. The equation set is obtained from the complete Navier- 
Stokes equations by retaining the viscous diffusion terms normal to the solid surfaces in every coordinate 
direction. For a body-fitted coordinate system (£, //. (j) fixed in time, these equations can be written in the 
conservative form as: 
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where U represents the conserved variable vector. The vectors F, G, H, and F v , G v , H v represent the 
convective and diffusive fluxes in the three transformed coordinate directions, respectively. In Eqn. (1), J 
represents the cell-volume or the Jacobian of coordinate transformation. A multigrid based, general puipose 
multi-block structured grid approach is used for the solution of the governing equations. In particular, 
the TLNS3D flow code is used in this study to solve Eqn. (1). Several references exist describing the 
discretization and algorithmic details of TLNS3D code. We include only a brief summary of the general 
features, and refer to the work of Vatsa and co-workers [1,2] for further details regarding the TLNS3D code. 


Numerical Discretization 

The spatial terms in Eqn. (1) arc discretized using a standard cell-centered finite-volume scheme. The con- 
vection terms arc discretized using second-order central differences with scalar/matrix artificial dissipation 
(second- and fourth- difference dissipation) added to suppress the odd-even decoupling and oscillations in 
the vicinity of shockwaves and stagnation points [3, 4]. The viscous terms are discretized with second-order 
accurate central difference formulas [1], For the present computations, the Spalart-Allmaras (SA) model [5] 
and the Menter’s SST model [6] arc used for simulating turbulent flows. 

For temporal discretization, the convective and diffusive terms arc grouped together, and Eqn. (1) can 
be rewritten as: 


— = — C(U) + D p ( U) + D a { U) (2) 

where C(U), D p ( U), and D a ( U ) arc the convection, physical diffusion, and artificial diffusion terms, 
respectively. The cell-volume or the Jacobian of coordinate transformation being included in these terms. 
The time -derivative term can be approximated to any desired order of accuracy by the Taylor series 
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Eqn. (3) represents a generalized backward difference scheme (BDF) in time, where the order of accuracy is 
determined by the choice of coefficients a 0 , a u a- 2 ... etc. For example, if ao = 1.5, a\ = —2 and a 2 = .5, 
it results in a second-order accurate scheme (BDF2) in time, which is the primary scheme chosen for this 
work due to its robustness and stability properties [7]. 

In Eqn. (3), the superscript n represents the time level at which the solution is known, and n+1 refers 
to the next time level to which the solution will be advanced. Similarly, n- 1 refers to the solution at one 
time level before the current solution. Regrouping the time-dependent terms and the original steady-state 
operator leads to the equation: 
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where E(\J n,n ~ 1 ’-) depends only at solution vector at time levels n and prior, and S represents the steady 
state operator or the right hand side of Eqn.(2). By adding a pseudo-time term, we can rewrite the above 
equation as: 
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Solution Algorithm 

The algorithm used here for solving unsteady flows relies heavily on the steady-state algorithm available 
in the TFNS3D code [1], The basic algorithm consists of a five-stage Runge-Kutta time-stepping scheme 
for advancing the solution in pseudo-time, until the solution converges to a steady-state. Efficiency of this 
algorithm is enhanced through the use of local time-stepping, residual smoothing and multigrid techniques 
developed for solving steady-state equations. In order to solve the time-dependent Navier-Stokes equations 
(Eqn. 5), we added another iteration loop in physical time outside the pseudo-time iteration loop in TLNS3D. 
For fixed values of ,S'(U" ), E(\J n,n ~ 1 ’"), we iterate on U" + 1 using the standard multigrid procedure of 
TLNS3D developed for steady-state, until desired level of convergence is achieved. This strategy, originally 
proposed by Jameson [8] for Euler equations and adapted for the TLNS3D viscous flow solver by Melson 
et. al [7] is popularly known as the dual time-stepping scheme for solving unsteady flows. The process is 
repeated until the desired number of time-steps have been completed. The details of the TLNS3D flow code 
for solving unsteady flows arc available in Refs. [7] and [9]. 


Boundary Conditions 

Except for the moving diaphragm, standard boundary conditions, such as the no-slip, no injection, fixed 
wall temperature or adiabatic wall, far- field and extrapolation conditions are used as appropriate for the 
various boundaries. The most accurate procedure to simulate the moving diaphragm would require moving 
grid capability. For simplicity, we chose to simulate this type of boundary condition by a periodic tran- 
spiration velocity condition. The frequency of the transpiration velocity at the diaphragm surface in the 
numerical simulation corresponded to the frequency of the oscillating diaphragm, and the peak velocity at 
the diaphragm surface was obtained from numerical iteration so as to match the experimentally measured 
peak velocity of the synthetic jet emanating from the slot exit. The pressure at the moving diaphragm is also 
required for closure. However, in the absence of unsteady pressure data from the experiment, we imposed 
a zero pressure gradient at the diaphragm boundary. We also tested the pressure gradient boundary condi- 
tion obtained from one-dimensional normal momentum equation [10], which had very little impact on the 
solutions. Due to the simplicity and robustness, we selected the zero pressure gradient boundary condition 
at the diaphragm surface. 
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Implementation and Case Specific Details 


Although the actuator geometry is highly three-dimensional, the outer flow field is nominally two-dimensional 
because of the high aspect ratio of the rectangular slot. For the present study, this configuration is modeled 
as a two-dimensional problem. A multi-block structured grid available at the CFDVAL2004 website is used 
as a baseline grid. 

The periodic motion of the diaphragm is simulated by specifying a sinusoidal velocity at the diaphragm 
surface with a frequency of 450 Hz., corresponding to the experimental setup. The amplitude is chosen so 
that the maximum Mach number at the jet exit is approximately 0.1, to replicate the experimental conditions. 
At the solid walls zero slip, zero injection, adiabatic temperature and zero pressure gradient conditions are 
imposed. In the external region, symmetry conditions arc imposed on the side (vertical) boundaries and 
far-field conditions arc imposed on the top boundary. A nominal free-stream Mach number of 0.001 is 
imposed in the free stream to simulate incompressible flow conditions in the TLNS3D code, which solves 
compressible flow equations. 

The code was run in unsteady (URANS) mode until the periodicity was established. The time-mean 
quantities were obtained by running the code for at least another 15 periods and averaging the flow quantities 
over these periods. The phase-locked average of flow quantities were assumed to be coincident with their 
values during the last full time period. 


Results/Discussion 

Only a brief set of results is included here, detailed results will be available in the workshop proceedings. 
Majority of the solutions were obtained with SA turbulence model on the baseline grid with a time-step of 
At = T per i 0 d/72, where represents the physical time for a complete period. In addition, results 

were generated for the following conditions: (a) coarser grid (eg) obtained by eliminating every other point 
from the baseline grid; (b) finer grid (fg) generated by increasing resolution in the normal direction by 50%; 
(c) lower time-step (lower dt), where At = T perto d/ 144 ; and (d) Menter’s SST turbulence model. The 
results for the vertical velocity near the slot exit center line (x=0., y=0.1 mm) from these computations are 
compared in Fig. 1, which shows very little sensitivity to turbulence model change, and spatial and temporal 
resolutions. 

Next we show the results for time-averaged velocity along the jet center line in Fig. 2. In the proximity 
of the slot exit (y < 2 mm ), all 5 sets of results are indistinguishable, and are in good agreement with the 
hotwire data. However, further away from the slot, the coarse grid (eg) results show significant deviation 
from the baseline results. The effect of turbulence model type (SA or SST) arc also more pronounced in the 
outer region. Similar trends arc observed in Fig. 3 for the jet width based on time-averaged solutions. 

Based on these results and more detailed comparisons (not shown here) of the computational results 
for time-averaged and phase-locked average velocities, the effect of lowering the time-step and refining the 
normal mesh were found to be very small. Noticeable differences were observed for the coarse grid (eg) 
results in the outer regions. The difference between the SA and SST turbulence model results become more 
pronounced in regions further away from the slot exit. 
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Figure 1: Time-variation of V- velocity at slot exit, x=0., y=0. 1 mm 
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Figure 2: Center-line average velocity comparisons 
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Figure 3: Jet-width comparisons 
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